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Abstract 

We investigate analytically and numerically the conditions for the Turing instability to occur in 
a one-dimensional chain of nonlinear oscillators coupled non-locally in such a way that the coupling 
strength decreases with the spatial distance as a power-law. A range parameter makes possible 
to cover the two limiting cases of local (nearest-neighbor) and a global (all-to-all) couplings. We 
consider an example from a non-linear auto-catalytic reaction-diffusion model. 
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I. INTRODUCTION 



Given the diversity of spatial patterns one observe in Nature, it is always puzzling to 
wonder how this diversity arises from spatially uniform situations. This question was pursued 

n 

by A. Turing in his seminal 1950's paper on the mechanism of morphogenesis in which he 
considered how can a spatially homogeneous state lose stability, giving origin to a spatially 
non-homogeneous pattern, the so-called Turing instability. Turing analyzed a reaction- 
diffusion system comprised of spatially coupled ordinary differential equations that modelled 
discrete biological cells. 

Turing put forward two basic assumptions: (i) the need of at least two chemical substances 
(morphogenes), an activator and an inhibitor; (ii) diffusion plays a destabilizing role in the 
interacting chemical substances. The latter point is somewhat counterintuitive, for diffusion 
usually smooths out spatial structures in linear systems and thus would be expected to 
exert a stabilizing influence instead. Moreover, Turing found that the instability caused 
by diffusion leads to the growth of spatial structure at a particular wave length. This 
mechanism has explained successfully a wealth of pattern formation phenomena in Nature, 
like the segmentation patterns in the developing fly embryo, the periodic array of tentacles 
around the mouth of the Hydra, zebra stripes sea shell striations [3], to name just a few. 
In physical chemistry Turing instability has been found in many continuously fed stirred 
tank reactors, like the Chlorine Dioxide-Iodine-Malonic Acid reaction 0, 5|. 

From the mathematical point of view, Turing instability arises when an otherwise stable 
homogeneous state becomes unstable as a result of diffusion, in the sense that any small 
perturbation will develop a spatially nonuniform state. Linear theory can predict the wave 
length of the unstable mode which grows 6|. This growth, however, is limited by the 
saturating effect of nonlinear terms and a spatially non-homogeneous pattern is then formed 
result of Turing instability. 

Turing instability is usually investigated in the context of local couplings, in which each 
cell interacts with its nearest-neighbors only. This is modelled by Pick's law, in which the 
diffusive flux of morphogenes (activator and inhibitor) is proportional to their local concen- 
tration, being responsible for the second spatial derivatives in the mathematical models of 
such reaction-diffusive systems. However, there are many situations in which this picture 
is over-simplified, like in chemical couplings between biological cells, generating non-local 



interactions which cannot be described by a diffusive couphng. Non-local couplings can be 
implemented in macroscopic density equations in which the diffusion coefficient depends on 
a weighted spatial average of the density 3] • 

Such nonlocal couplings take into account not only the nearest neighbors of a given cell, 
but also the other cells of the assembly. An extreme form of this arises when each cell 
couples with the mean field produced by all the other cells, the so-called global (all-to-all) 
coupling, which is used in models of neural networks js], coupled Josephson junctions jo]. 



and it is the basis of the paradigmatic Kuramoto model 10|. It is possible to interpolate 



between nearest-neighbor and all-to-all coupling by considering an interaction term whose 



strength decreases with the physical distance between cells as a power-law 



This model 



has been used to theoretical and numerical studies of various spatio-temporal phenomena 



like synchronization, shadowing and cluster formation 12|. 



The influence of non-local couplings on the formation of Turing patterns has been studied 
in many systems. The conditions for the occurrence of Turing instability were investigated 
for logistic growth processes using various non-local coupling terms [l3|. In other recent 
study, a three-variable Oregonator model of a light-sensitive Belousov-Zhabotinsky reaction 
was studied, in which the nonlocal coupling was externally imposed by an optical feedback 
loop [l4|. In the latter investigation, it was found that long-range inhibition leads to Turing 
instability, whereas long-range activation induces wave-like patterns. 

This work aims to investigate the presence of Turing instabilities in oscillator chains with 
a non-local coupling of the power-law form, so as to have results which hold both for local 
and global cases. We derive analytically the conditions for a Turing instability to occur in 
a one-dimensional lattice of two-dimensional systems (representing the activator-inhibitor 
pair) undergoing linear dynamics. Then we study a nonlinear reaction-diffusion system 
proposed by Meinhardt and Gierer as a model for pattern formation related to skin 
pigmentation. 

The rest of this paper is organized as follows. In Section II we introduce the oscillator 
model to be studied, as well as the type of coupling between the units. The analytical 
conditions for a Turing instability to occur are derived in Section III using linear stability 
theory. The effect of nonlinear terms and the ensuing pattern formation is discussed in 
Section IV through the Meinhardt- Gierer model. Our conclusions are left to the last Section. 



II. OSCILLATOR CHAIN WITH NON-LOCAL COUPLING 



Let us consider an oscillator one- dimensional chain with units, each of them being a 
nonlinear react ion- diffusion system, whose dynamical state is described by the concentrations 
of the activator Xkit) and inhibitor yk(t), satisfying the following equations {k = 1,2, ■ ■ ■ N): 

N' 



Xk = X(a;fc, yk) - D^Xk + ^ {xk-r + Xk+r) , (1) 

D ^' 1 

Vk = Y(Xfc, Vk) - DyVk + — ^ V — iVk-r + Vk+r) , (2) 

r=l 



where X and Y stand for the vector field corresponding to the uncoupled oscillators. We 
use a non-local form of coupling for which the interaction strength decreases with the lattice 
distance in a power-law fashion, where and Dy are positive coupling constants, repre- 
senting the different diffusion coefficients of the chemical species, a is a positive real number, 
and 

TV' ^ 
r=l 

is a normalization factor, with A^' = (A^ — l)/2, supposing that is an odd number. In 
the oscillator chain given by Eqs. ([I])-®, the summation terms are weighted averages of 
discretized second spatial derivatives, the common normalization factor K{a) being the sum 
of the corresponding weights. 

In the limit a — )■ oo, only the r = 1 term will survive in the summation terms, resulting 
in K — 2 and the Laplacian or diffusive coupling 

Xk = X(xfc, Vk) + ^ (xk-i - 2xk + Xk+i) , (4) 

yk = Y{xk, Vk) + ^ iVk-i - 2yk + yk+i) , (5) 

in which only the nearest-neighbors of a given site contribute to the coupling term. The 
other limiting case, a = 0, is such that k(0) = A^ — 1 and the oscillator chain becomes 
globally coupled 

Xk = ^{xk,yk) + Dx {-Xk + j^j—j ^ Xr\, (6) 

\ r=l,ry^k / 



Vk = Y{xk, Vk) + Dy I -yk + ) ' (7) 

r=l,rj^k J 



where each oscillator interacts with the mean value of all lattice sites, irrespective of their 
positions ( "mean- field" model). Hence, the coupling term in Eqs. (II])-© may be regarded 
as an interpolating form between these limiting cases, and will be referred to as a power-law 
coupling. 



III. LINEAR STABILITY 



Let us consider the equilibrium values of the activator and inhibitor given, respectively, 
by Xq and Yq. The deviations from these values (also denoted by x and y, for notational 
simplicity), in lowest order, obey a linearized model, for which the activator-inhibitor dy- 
namics in each uncoupled oscillator (labelled by k, as before) is governed by the affine vector 
field 

^{xk,yk) = axk + byk, (8) 
Y(xk,yk) = cxk + dyk, (9) 

where the constants a to d are the elements of the Jacobian matrix of the nonlinear vector 
field, evaluated at (Xq, Yq). The equilibrium solution of the linearized system is (x* = 0,y* = 
0), which is linearly stable if 

q = ad — be > 0, and a + d < 0. (10) 

The convergence to the equilibrium will be monotonic if (a — rf)^ + 46c > and oscillatory 
if this factor is negative. This set of conditions shall be supposed from now on. 

In order to treat the coupled system we can use a discrete Fourier transform of each 
dynamical variable, in the following form 



Xk{t) = Y,Us,t)exp [-^j , (11) 



s=0 
N-l 



Vkit) = J2vs{s,t)exp (^^) , (12) 

where and rjs are time-dependent Fourier mode amplitudes. On substituting (fTTl) - (fT2l) in 
the linearized version of dl])-® we find that the mode amplitudes satisfy uncoupled linear 
differential equations: 

is = [a- 2D^cr„(s, N)]^s + br]s = a^^s + br]s, (13) 
fi^ = [d- 2Dyaa{s, N)]^s + ct], = c^^ + d^r]s, (14) 
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FIG. 1. Dependence of the function defined by Eq. (jl7p on the argument s/N, for (a) small and 
(b) large values of the parameter a. 



where we define the following auxiliary quantities: 

da = d — 2Dy(Ta{s, N), 

N' 



2 K,(a) ^ 



cos 



2T\sr 
N 



(15) 
(16) 

(17) 



is a function whose values are such that < crQ,(s,A^) < amax{<^, N) < 1 [Fig. [T]. For the 
nearest-neighbor case (a — )■ oo) it turns out that 



croo(s, A^) = sin^ 



2 

n) 



(18) 



for which (Tmax = 1 reaches its largest value [Fig. lU^b)]. On the other hand, in the mean-field 
coupling (a = 0) this function reads 



I f N 



2\N -I 



(19) 



such that, for large it tends to 1/2 [Fig. [U^a)]. 

The Turing instability occurs when a spatially homogeneous pattern becomes inhomo- 
geneous by virtue of the coupling effect which acts as a perturbation on each oscillator, 



driving it out of its equilibrium state. The equilibrium point of the linearized mode am- 
plitude equations, {^s,Vs) = (0,0), corresponding to equilibrium concentrations of activator 
and inhibitor species for each oscillator. Hence the onset of a Turing instability happens 
whenever this equilibrium becomes an unstable saddle point, which occurs if 

= a^d„ -bc = AD^Dya^ - 2a{aDy - dD^) + g < 0. (20) 

On defining the auxiliary variables 

P.^,^, pi) 

the inequality (120|) is satisfied provided a_ < cr < cr+, where 

cy± ^ ^4 -. (22) 

Since the function o has the upper bound Omaxis^-,^^ < 1 we must have < (t_ < 
<ymax{p-^ N), such that < P — a/P^ — 4Q < 4amax{(^, N), which gives the following condi- 
tions for the appearance of Turing instability in the system 

g > 0, (23) 



P>2^Q, if 0<P<4or„,,,(a,iV), (24) 
P>n T— r77 + 2a^ax(«,iV), if P>Aa^ax{a,N). (25) 

Hence the stability of the spatially homogeneous state in this system depends on the 
values taken on by <Jmax- We show in Fig. [2] how it depends on the range parameter a by 
direct numerical evaluation of the function extremum (solid curve) and using an analytical 
approximation (dashed curve). The latter was obtained as follows: first we can show from 
f|T7|) that this function is always symmetric with respect to s = N/2: 

a^{s = N/2 + L,N) = a^{s = N/2 -L,N), 

for all L values such that < s < — 1. 

From the same token we show that aa{s,N) has a local maximum at s = N/2 (for 
N' = {N-l)/2 odd) 

dermis, N) 



ds 



s=N/2 



< 0. 

s=N/2 
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FIG. 2. Dependence of the maximum value of aa{s,N) on the parameter a. The full curve 
represents a numerical evaluation of the summation, whereas the dashed curve results from an 
analytical approximation. 

such that 

(Tmax{a,N) = a^{s = N/2,N). 
Finally, if is large enough that s/N can be reasonably well approximated by a con- 



tinuous variable, we can use the Euler-MacLaurin formula 



16| to replace the summation in 



f fT7|) by an integral, which yields an analytical approximation for a„ 



<Jr. 



(«,iV) 



1 - {N') 



a 



1 



la— 1/ AT' 







I fl- (Nf 


2 




2 I a 


- 1 




{N' - 


l)"" + 2-"" 


)■ 




2 





N'-" + 1 



(26) 



a — 1 

which holds for odd values of A^'. We see that the approximation always underestimate the 
actual value of amax{(^, N), but the difference turns not to be large, specially at small a. 
On the other hand, for large a, while the analytical approximation gives worse results, we 
already know that this value must approach the unity, since this is the nearest-neighbor 
coupling case. Notice that, in this work, the large- iV limit has to be taken only as an 
approximation for obtention of closed- form analytical expressions. In fact, for — oo the 



power-law coupling ([I])-(l2]) is not normalizable when < a < 1 17l |. 

In the case of local (nearest-neighbor) coupling {a — t- oo) we have large values of a and 
<^max tends to the unity, such that the conditions (|23|) - (|25|1 a Turing instability to occur 
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reduce to 



g > 0, (27) 



P > 2y/Q, if < P < 4, (28) 
P > I + 2, if P > 4, (29) 



which agrees with the results of Ref. |6| up to a (nonessential) factor of 2 due to a slightly 
different definition of the local coupling prescription. 

Since he uncoupled oscillators are supposed to have a stable equilibrium concentration, the 
local conditions (fTOl) apply, such that ( 1271) always holds. The remaining conditions depend 
on the values of the diffusion coefficients. It is instructive to consider the case where both 
diffusion coefficients are equal: = Dy = D. In this case it turns out that P = {a + d)/D 
is negative, and there cannot occur a Turing instability. In fact, it is a well-known fact that 
pattern formation in a chemical system will not occur unless the diffusion coefficients of the 
activator and inhibitor differ substantially jsl. 

The range of unstable modes, in the local coupling case, is given by 

min(cT+, 1) > a > a_, (30) 

where a± are given by ( l22i) . Notice that we have ct+ > 1 provided P < {Q/2) + 2. Substitut- 
ing (|T8l) and solving the resulting inequalities there result two intervals of unstable modes 
symmetrically located with respect to s/N = 1/2: 

sin""*^ ^/aZ s sin~^ ^/al 



< ^ < ' 

1 ^ < T7 < 1 — , (31 

In this way, a larger number of modes become unstable if 0"+ > 1. 

The global (all-to-all) case (a = 0) presents interesting features, for the function cro(s, A^) 
takes on a constant value ctq given by (fT9l) . The conditions ( I23|) - (p5l) hence reduce to a single 
inequality 

P > ^ + 2ao, (32) 
2 (Jo 

in such a way that, once this condition is fulfilled, all modes become simultaneously unstable. 
Since, for large A^, we have ctq ~ 1/2, the relation above can be written in a simple form P > 
Q+l. It is interesting to note that, if the diffusion coefficients are equal, this condition cannot 




FIG. 3. (a) Parameter plane showing the hnes of marginal stability for the case of = 101 coupled 
oscillators and arbitrary a. (b) Magnification of a region of (a) of large Q-values 

be satisfied in general, as in the local case. Thus there is a range- independent influence of 
the dissimilarity of the diffusion coefficients on the occurrence of a Turing instability. 

Let us consider a specific example of the intermediate range case, for which a = 1 and 
we take a chain with = 101 oscillators. Then (IT7|) gives (Jmax ~ 0.577 for this case. The 
approximation expressed by Eq. (!26|) . which works well also for A^' even, would give a slightly 
smaller value of 0.568 for cTmax- On substituting the value of (Tmax in Eqs. f l23|) -f l25|) we obtain 
conditions for the existence of a Turing instability which can be better represented in the 
parameter plane depicted (for a arbitrary) in Fig. [3]^a), where we draw the lines where we 
have marginal stability as a function of parameters P and Q. Since, for a = 1, the function 
a has a broad flat top, failing on being a constant (like in the global case) only for the 
vicinity of the extremities s/N = and s/N = 1, we have roughly the same behavior in 
that most modes become unstable simultaneously in the intermediate range case. 

A useful information is the stable area Agt, or the area in the parameter plane for which 
we have stable equilibria. In Fig. El^b) we show a magnification of a region of large Q- 
values in the parameter plane for different values of a. It is apparent that the stable area 
decreases with increasing a, i.e. the relative number of parameter values yielding stable 
equilibria is larger as the coupling becomes global, and smaller as the coupling becomes 
local. By integrating out the stability curves in Fig. ^a), from Q = to an arbitrary Qo, 
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it is actually possible to derive an analytical expression for the stable area (normalized by 
(5o)j which reads 



^ ^ ^ 1 ^ ^max(2Qo-4) ^ 'i2oVL _ ^^^^ 



giving, for large Qo? ~ l/(4crma2,). In very long oscillator chains, where is large, with 
local and global couplings a max takes on values equal to 1 and 1/2, respectively. Hence the 
stable area in the local case is roughly half of the corresponding area in the global case. We 
thus conclude that Turing instability would be statistically more common for local couplings 
than for global ones. 

IV. NONLINEAR REACTION-DIFFUSION MODEL 

In reaction-diffusion systems of practical interest, there are nonlinear terms which exert 
saturation on the non-homogeneous spatial modes excited by the Turing instability, giving 
rise to pattern formation. The basic mechanism of this saturation combines a short-ranged 
auto-catalytic activation, which causes self-enhancement after Turing instability sets in, with 
a long-ranged inhibitory effect. The latter can be caused by either a rapidly spreading and 
long-ranging or a result from a depletion of material required for the self-enhancement that 



is obtained from the surrounding region 15 1 



One paradigmatic example of this broad class of react ion- diffusion 



Meinhardt-Gierer model which, in the language of Eqs. ([I])-®, reads 



15 



systems is the 



id 



X(xfc, yk) = px— - PxXk, (34) 

Uk 

Y(xfc, yk) = pyxl - PyVk, (35) 

where pa^h and pa,h are positive parameters characterizing the local reaction rates. 

The term proportional to p^ comes from the assumption that the activator is an auto- 
catalytic chemical species, its denominator representing the inhibitor effect of the other 
species. In this case p^ turns to be a small activator-independent production rate of the 
activator and is actually necessary to start the activator auto-catalysis at very low activator 
concentration, like in the case of regeneration. Likewise, py stands for a baseline production 
rate of the inhibitor 1^. The quantities p^ and Py play the role of degradation of the 
activator and inhibitor species, respectively, such that the number of molecules of each 

11 



species that decay per time unit is p roportional to the corresponding decay rate and to the 
number of molecules themselves 19 1. 



A. Stability analysis 

We first determine the equilibrium points of the uncoupled vector field (1M|) - (!35|) . One of 
them is obviously x* = y* = 0, whereas there is a non-trivial equilibrium with activator and 
inhibitor concentrations given, respectively, by 

X* = y* = (36) 

Pypx PyPx 

The linearized dynamics around this equilibrium, for the uncoupled case, is obtained from 
the Jacobian matrix of ( !34|l - (!35|l . with elements evaluated at the point {x*,y*), giving the 
matrix elements 

a = Px, b = -^, (37) 

Px 

c = d = -py, (38) 

Px 

If the equilibrium {x*,y*) is to be stable, the trace and determinant of this Jacobian must 
be, respectively, negative and positive, yielding the following conditions 

Px < Py, PxPy > (39) 

Since all coefficients are supposed positive, this implies that the decay rate of the inhibitor 
must be greater for the activator, which is a reasonable assumption for having a stable 
equilibrium concentration of both species. The approach to the equilibrium is determined 
by the nature (real or complex) of the eigenvalues of the Jacobian. The results of this simple 
analysis are summarized in Fig. |H By the same token, the equilibrium at the origin can be 
shown to be always unstable and thus it is not interesting when discussing Turing instability 
due to diffusion. 

Now we turn into the coupled system ([I])-®, for which there is an spatially homogeneous 
equihbrium state given by 

Xk=x*, yk = y\ (A; = l,2,---iV) (40) 

This state defines a two-dimensional invariant subspace embedded into the full 2A'^-dimensional 
phase space of the coupled oscillator chain. The subspace is invariant for any trajectory 
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FIG. 4. Stability of the uncoupled equilibrium points of the Meinhardt-Gierer equations in the 
parameter plane [i^ versus fiy. 

which originates from an initial condition belonging to this subspace is bound to lie within 
it for further times. Geometrically speaking, the possible spatially non-homogeneous states 
one may find represent trajectories outside this invariant subspace. 

In this geometrical language, we are interested to investigate the stability of the invariant 
homogeneous subspace along the 2{N — 1) remaining transversal directions. We can use the 
results of Section II for obtaining analytical conditions for the stability of the spatially ho- 
mogeneous state. Substituting (!37j) - (!38l) into (I2T]) there results that the stability conditions 
(El-fESD read 



(41) 



II^Dy - HyD^ > fl^flyD^D. 



2a max 



yi 



if < yUa;-Dy - /iy-Dj; < 4:a.maxDxDy, (42) 

if ^.^Dy - fiyDr, > 4(7 maX D X D y , (43) 



Since the uncoupled oscillators are supposed to be locally stable, so that Eq. fl5Ul) holds, the 
condition f HT]) is always fulfilled, since the diffusion coefficients are positive-definite. 

Let us consider the limit cases of this coupling. For nearest-neighbor case (large a) we 
have (Tmax = 1. If the diffusion coefficients were equal, there would be no Turing instability, 
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0.01 




FIG. 5. Parameter plane of the activator and inhibitor diffusion coefficients for = 0.01, Hy = 
0.02, px = 0.01, and py = 0.02. 

as shown before. For the all-to-all coupling (a = 0) the conditions (H2l)-(H3ll reduce to the 
single inequality 

HxDy - ^yDx > fixfJ-y + DxDy (44) 

The threshold for Turing instability is depicted in the parameter plane versus Dy (for 
fixed values of the remaining parameters) [Fig. |5] . 



B. Pattern formation 



In the numerical simulations to be presented, we use a predictor-corrector routine 



(LSODA) based on the Adams method 20|], and the same values for the parameters 



of the Meinhardt-Gierer model as those chosen in Ref. [2l|: Hx = 0.01, fiy = 0.02, 
Px = 0.01, and py = 0.02, for which the equilibrium point is a stable focus with coordinates 
{x*, y*) = (1.0, 1.0). We also fix the diffusion coefficient of the inhibitor at Dy = 0.2 and use 
the activator coefficient Dx as the tunable parameter for each value chosen for the effective 
range a. We use, as an initial condition, an almost homogeneous spatial profile with a 
tiny bump which is intended to be a seed for some unstable mode to grow, after a Turing 
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instability occurs. 

The local case (large a) is illustrated by Fig. [6l where we depict space-time diagrams 
(activator concentration x versus discrete space k and time t) for different values of the 
activator diffusion coefficient. For relatively large [Fig. [H^a)] we do not observe pattern 
formation, since we are ahead from the instability threshold predicted by the linear theory. 
The oscillators go together to the homogeneous equilibrium state x* = 1 instead. The 
instability threshold occurs for ~ 0.016 and we accordingly observe the emergence of a 
spatially non-homogeneous pattern just after the Turing instability [Fig. [6]^b)]. 

This pattern is roughly a sinusoidal spatial profile with 12 maxima, corresponding to a 
wave-length of A ~ 8.42 and thus to a wave number 

From the linear stability analysis, the parameter set used to generate Fig. Et^b) corresponds 
to cr_ = 0.09 and a+ = 0.17. Using Eq. ( I3T1) . the ranges of unstable modes are [0.097, 0.135] 
and [0.865,0.903]. The pattern observed has wave number roughly in the middle of the 
interval of unstable modes, i.e. it is the most unstable mode from the linear approximation. 

As we decrease substantially the activator diffusion coefficient, we observe another spatial 
pattern [Fig. EJ^c)] with 14 maxima and thus a smaller wavelength than the previous example, 
namely 

^ ^ 0.139, 



N 7.21 

which belongs again to one of the intervals of linearly unstable modes: [0.076, 0.392] and 
[0.608,0.924]. It is no longer the most unstable mode, however. A larger wave number 
such as this, is actually a general feature of spatio-temporal systems when we decrease the 
diffusion coefficient: modes with small wave number are more easily damped when diffusion 
increases, as a general rule. 

The corresponding space-time plots for the global (all-to-all) coupling case are shown in 
Fig. Wi^) and (b) for parameters before and after the occurrence of a Turing instability, 
respectively. According to the linear theory [see Fig. \5\ it occurs for ^ 0.008. The 
pattern formed after the Turing instability is more complex than in the previous example 
and cannot be assigned to a single unstable mode, but rather to a entire spectrum of unstable 
modes excited by the Turing instability, and which became stabilized by the saturating effect 
of the nonlinear terms of the model, a feature obviously not present in the linear stability 
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FIG. 7. (color online) Space-time diagrams of the activator concentration for the all-to-all (global) 
coupling case, = 0.2, and (a) Dy = 0.009; (b) Dy = 0.005. 




FIG. 8. Space-time diagrams of the activator concentration for the intermediate range coupling of 
a = 1.0, = 0.2, and (a) Dy = 0.020; (b) Dy = 0.005 

analysis of the previous Section. The space time plots for the intermediate range coupling 
characterized by a = 1.0 are depicted in Fig. M^a) and (b), before and after the threshold 
of a Turing instability, respectively. 
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V. CONCLUSIONS 



Turing instability is a general mechanism of pattern formation, from which a spatially 
homogeneous pattern becomes unstable and develops a spatial structure after saturation 
that comes from nonlinear effects. In biological scenarios there may be the need for nonlocal 
couplings, or interactions among distant cells, not only the nearest neighbors, as in usual 
studies of spatial diffusion. As far as we know, this is the first analytical study of the 
occurrence of Turing instability in coupled oscillator chains with nonlocal interactions. We 
presented a linear stability analysis for the Turing instability when the coupling among 
sites is non-local, i.e. it takes into account not only the nearest-neighbors but also all the 
other oscillators in a one-dimensional chain, the coupling strength decreasing with the lattice 
distance as a power law. 

We shown, using a specific nonlinear model, that the linear stability analysis is useful to 
characterize the Turing instability when nonlocal couplings are considered. In particular, we 
derived analytical conditions for the presence of Turing instability for an arbitrary effective 
coupling range. Our results agree, in the limit of short range, with results obtained for 
nearest-neighbor couplings, and can be also applied to situations with global (all-to-all) 
couplings. We shown that local couplings have a comparatively shorter stable area in the 
parameter plane than global couplings. In other words, as the effective coupling range 
increases, the more stable are the spatial configurations and, statistically speaking, less 
probable would be the occurrence of a Turing instability. 

This fact can be qualitatively understood by noting that global couplings spread infor- 
mation among cells more rapidly than locally coupled ones, for which diffusion happens at 
a slower rate. Then globally coupled cells are less likely to present Turing instability than 
locally coupled ones. On the other hand, other dynamical collective phenomena due to in- 



teractions, like frequency synchronization [22], short-term memory formation 23|, bursting 



synchronization 



24j |. are more likely to occur in the globally than locally coupled chains. 



The question of pattern formation, however, is more complex, since the patterns are 
ultimately determined by nonlinear features of the model not predictable by linear stability 
analysis. Even so, in some cases it is possible to relate single-mode patterns to the modes 
which become unstable after a Turing instability. 
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